The noise propagator in an optical 
system using EDFAs and its effect on 
system performance: accurate evaluation 
based on linear perturbation 

Zhongxi Zhang*, Liang Chen, Xiaoyi Bao 

Department of Physics, University of Ottawa, 1 50 Louis Pasteur, Ottawa, K1N6N5, Canada 

* zzhan3 @ uottawa. ca 

Abstract: For optical systems amplified by EDFAs, we present a new 
approach to obtain the noise propagator from the accurate linearized noise 
equation (LNE) proposed by Holzlohner et al (Tj. With the help of fourth- 
order Runge-Kutta in the interaction picture (RK4IP) method, we solve both 
the accurate LNE and the nonlinear Schrodinger equation (NLSE) for the 
noise-free signal simultaneously. The noise propagator matrix thus obtained 
does not need the jitter separation required by the covariance matrix method 
(CMM) U. Based on this new approach, the calculated BER in a DPSK 
system agrees well with the simulation using CMM. By introducing a 
deterministic phase shift induced by Gordon-Mollenauer noise effect, our 
numerical results also agree well with the experimental data of multi-span 
DPSK systems. 
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1. Introduction 

The amplified spontaneous emission (ASE) noise generated from optical amplifiers, e.g., Er- 
bium Doped Fiber Amplifiers (EDFAs), is one of the basic reasons for the bit-error-rate (BER) 
in an optical fiber communication (OFC) system. For an OFC system with non-negligible Kerr 
nonlinearity, the ASE impact evaluation is complicated because of the nonlinear interaction 
between signal and ASE noise. In the case of signal field being much stronger than the noise 
field, the noise-noise beating is relatively small so that the noise field in the fiber can be approx- 
imately described by an equation called linearized noise equation (LNE), which was proposed 
separately in Ref. Q] and Refs. CJ-O. 

Noise propagator is a matrix introduced to describe the noise propagation in the fiber. Within 
the framework of linear perturbation, it is independent of the specific noise realizations and can 
be used to evaluate the noise impact on the moment-generating function (MGF) of the filtered 
photoelectric current and then on the BER at the receiver IT] [5] El El - 

The popular form of LNE (2|-[8| is based on the continuous wave (CW) assumption, i.e., 
the noise-free signal in the LNE is artificially simplified as a CW wave. As a result, a semi- 
analytical form of noise propagator and the noise power spectral density (PSD), usually called 
parametric gain (PG), can be obtained. The drawback of this simplification is that the noise-free 
signal in this LNE has no chromatic dispersion (CD) effect. Thus the couplings between noise 
components (in frequency domain) cannot be taken into account. 

A LNE beyond CW, named accurate LNE in this work, was first proposed and discussed in 
Ref. UJ . Dynamically taking into account the local CD and Kerr nonlinearity along the fiber, 
this LNE provides accurate noise information, with its computational cost being much higher 
than the CW approach. As an example, given a noise-free signal obtained from NLSE, the com- 
putation required to update the accurate LNE has cubic complexity in the number of Fourier 
components [1]. To reduce the computational complexity, co variance matrix method (CMM) 
was proposed in Ref. JT), where the noise co variance matrix was obtained by processing large 



enough noise realizations. In Refs. Il9l [Toll , the computational cost of CMM was further re- 
duced by a deterministic approach using perturbation solution. Since the raw covariance matrix 
obtained either by noise realizations JT] or by deterministic approach |9][l0) may contain non- 
linear noise contribution, it is important to separate the nonlinearity-induced phase jitter and, in 
some cases, time jitter from the raw covariance matrix. Thus, the obtained pdfs of the receiver 
voltage agrees well with Monte Carlo simulation lfTll9l [TT1 . 

Based on linear perturbation, the noise covariance matrix can also be obtained by solving its 
ordinary differential equation (ODE) proposed by Q][l2). The covariance matrix obtained by 
solving such linear ODE does not need jitter separation, although this ODE is more complicated 
than the accurate LNE lfl2l . So far there is little comparison between the approaches of Ref. 
Ifl2l and Refs. [QIUCEO). 

In this work, we show that the noise propagator matrix can be obtained directly from the 
accurate LNE. Therefore, there is no nonlinearity-induced jitter in such noise propagator. To 
effectively reduce the computational complexity in updating the LNE, one can decompose the 
Kerr effect related matrix into a symmetric and an antisymmetric matrices [cf. the discussion 
after Eq. d36b in the Appendix]. Also, to solve the accurate LNE with large enough step size, 
the fourth-order Runge-Kutta in the interaction picture (RK4IP) method lfT3l[T4l is generalized 
and applied, as shown in Sec. [2] Based on Refs. fl~|[5]0|8), we evaluate the impacts of noise 
propagator on MGF and BER in Sec. [3] The accuracy of this new approach depends on how far 
the linearized noise deviates from the actual noise. To numerically verify this new approach, we 
calculate the BERs in a 20-span DPSK system discussed in Ref. Q, where the CMM results 
with nonlinear phase noise = 0.2n were given. Our BER calculations agree well with the 
CMM results. We also consider the experimental data of the multi-span DPSK systems in Ref. 
JT) ■ A new phase shift is introduced to fit the experimental data. This new deterministic shift is 
determined by the variance of the Gordon-Mollenauer (GM) phase noise lfl5ll . 

2. Noise propagator obtained from the accurate LNE 

Within the framework of linear perturbation, the noise propagator in an EDFA-based system is 
a fundamental matrix that determines the noise impacts on MGF and BER. In this section, we 
show that the noise propagator matrix in a fiber of length L can be obtained from the accurate 
LNE given by Eq. (I38t . For a multi-span OFC system, one needs to introduce an equivalent 
noise propagator which can be obtained by calculating the PG. 

2.1. Noise propagator in a fiber of length L 

The noise propagator in a fiber of length L can be obtained by extending the RK4IP in Refs. lfl3l 
[141 to the LNE ([38]). By introducing a = e^^a' and = e ~ L{ - z ^Ne L{ - z ^ = e' t( - z ~^ (v + 
(i)e L ( z ~ Z0 \ Eq. d38l in the interaction picture (IP) has the form 

f=N'a' (1) 
dz 

Taking zo = z„ + h/2 with h = z„+i —Z n and denoting a n = a(z„), a n+ \ = a(z n +l), Nj, = 
e f "'^ 2 N(zn)e- M / 2 , N' n+lI1 = N(z„ +h/2), and N' n+l = e- Lh l 2 N{z n +\)e lh l 2 , one can use RK4IP 
inUHl to to solve Eq. © with 
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which means the noise propagator for the fiber of length h = z n +\ 



^N{z n+i )e Lh l 2 {e Lh l 2 + hh)\a n , 



z„ can be calculated as 
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For the fiber of length L, the noise propagator has the form 

p n {L,0)=H(L,L-h L )---H{h u O). 



(2) 
(3) 

(4) 
(5) 



Note that the RK4IP used here is different from the RK4IP in Ref. 1 14], where what to be 
solved was the noise-free signal (a ID matrix), while here what we want is the noise propagator 
(a 2D matrix), whose computational complexity is 0(N^,), due to that each kj (i = 2,3,4) in 
Eq. d2J needs one general matrix multiplication. 

2.2. Equivalent noise propagator of a multi-span system 

As discussed in the Appendix, the ASE from an EDFA can be modeled as additive white Gaus- 
sian noise (AWGN) with its variance given by Eq. ( 1391 . Given the ASE injected at the input of 
a fiber and the noise propagator obtained from Eqs. (O, (0J, and (0, the noise PSD (or PG) at 
the output of a fiber of length L can be written as U |5] H [8] 



Pd = Pn (L,0)a 2 I P l(L 1 0) = o 2 Pl1 (L,0)pl(L,0), 

where / is a unit matrix and Eq. d39t has been used. In Eq. (0, p T n is the transpose of p„. 
For a K-span system consisting of (K + 1) EDFAs (cf. Fig. |4j), its PG has the form 



(6) 



PG = (Ga 2 + a 2 )P„(K, 0)Pl (K, 0) + a 2 £ P n (K, k)P* (K, k) 



k=\ 



P n (K,k) = p„ (KL, (K-l)Lj ■ ■■p n ({k+\)L,kL^ , (k = Q,l,---,K—l) 



(7) 



where af n = N{ n / (2Tq) and a 2 is given by Eqs. (l28l i and 
matrix PG is positive definite. It can be factorized as 



PG (7 Pn.eqPji^eqi 



In Eq. (O, the real symmetric 



(8) 



where the equivalent noise propagator P ne „ can be obtained either by using Cholesky decom- 
position or symmetric (square root) decomposition ifTBI . The later yields P n ,eq = Pneq- O ur 
calculation confirms that both decompositions result in the same BER. Note that the P n . eq ob- 
tained from Eq. ([8]) does not necessarily mean it can provide all the noise information. For 
example, the real unitary matrix associated with a phase shift in noise will lead to the same 



PG. This phase shift in noise will not affect the signal-signal beating and noise-noise beating. 
But it will affect the signal-noise beating. Our numerical results in subsection 15 .21 show that a 
deterministic phase shift in noise will significantly affect the calculated BER. Based on the ex- 
perimental data, we find that this phase shift is related with the variance of Gordon-Mollenauer 
phase noise. Here we name it the GM shift Aqm- 



3. BER calculation 

With the noise propagator matrix obtained from the accurate LNE, one can evaluate the BER 
in the OFC system by calculating the moment-generating function (MGF) of the electrically 
filtered current 7(f s ) using Karhunen-Loeve series expansion (KLSE). 

3.1. Filtered photoelectric current expressed using KLSE 

Given a linear optical system, based on the discussions in Ref. ifTTl and the notations introduced 
in Ref. fT8l , the filtered photoelectric current 7(f) can be expressed in the form of 7(f) = \(s" (t + 
T b )+n"(t + T b )\s°(t)+n"(t)) +c.c.]/2. Here s°(t) (n"(t)) represents the signal (noise) field at 
the input of the optical filter. Dirac bra (x\ is the conjugate transpose (or Hermitian transpose) 
of Dirac ket |x) [x = s°(t),n°(t),s°(t) +n°(t), etc.]. The Dirac ket differs from usual complex 
vector in that the z'th element of the later is just the z'th Fourier coefficient of the (signal or noise) 
field, while the z'th element of the former is the product of the z'th Fourier coefficient and its base 
function (cf. Eq. (17) in Ref. |fl8l ). According to Refs. |fT71[T8l , the filtered current in a linear 
system can be formally expressed as 7(f) = y ss +y n „ +y„ s with (/ = — L s , ■ ■ ■ L s ; m — ~M„ ■ ■ ■ M n ) 

y„(h) = [{s°(t s + T b )\R ss \s°(t s ))+c.c.]/2= (s°(t s )\R°\s°(t s )) 
Junih) - [(n (t s + T B )\R nn \n (t s )}+c.c.]/2=(N°\RZ\N°) = (Z\A\Z) 
y«(fi) = [{n°(t s + T b )\R m \s"(t s )) + (n"(t s )\R m \s°(t s + T h )) +c.c.}/2 

= l(N in \OlR°V(h)) +c.c] = [(Z\b D (t s )) +C.C.} (9) 

where A = U* Olfi'^O^U = diag{X u - ■ ■ ,X 2 m„+i}, \b D (t s )) = U^Ol n R^ s \s°{t s )), and 

(Rss)l'l = (^.w)/'/^/' 1 /, ( R nn)m'm = {Rnn)m'mD™ m , (R® s )ml = (Rns)mlD'ml 
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In Eqs. (t9T>-(fT0b. |Z) represents the decoupled Gaussian random variables with zero mean and 
real part and imaginary part variance of a 2 . The effects of the optical and electrical filters in 
the receiver are represented by matrices with their elements being (<9„„ ) mm i = o~ mm tH {jr), 



{O ss )w = Su'Hoim;) and (R SS ) IV = H T {$jfi), (7? ; 



H r {^),(R, 
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Due to the optical and electrical filters, signal (noise) components outside ±L S (±M„) can be 
neglected. Here IfTTl 



L s = r]NT b B a , M n = 7]B o T , T Q = n( 



(11) 



B () B,' 

For a nonlinear optical system, to get the noise propagator from the accurate LNE (1351) . one 
needs to separate complex numbers into their real and imaginary parts. Denoting the Re-Im 

Re{x} — Im{x} 
Im{x} Re{x} 

in Eq. ©J and introducing \n") = P n ,eq\ao) with \ao) being the AWGN from EDFA, 



form of a complex matrix x as x = 



[where x can be any complex matrix 
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it is easy to generalize the noise related currents in Eq. (O as |[Tll5ll7l[8l [T7l 

y nn {h)={a a \UAU T \a ) = {Z\R\Z) (A=U T Pl eq dZ n R° n d nn P nteq U=diag{l u - ■ • ,A 4M „ +2 }) 
y„(u)=^\0 T P^dl^6„\^h))=(^\B\nh)) = (Z\b(t s )) (B=U T Pl ec) dJ m R°d ss ), (13) 

where Pn,eq can be obtained from Eqs. (|2]i-(|5]l and C7}-®- 

3.2. BER obtained from MGF 

Averaging |Z)j (i =!,-■■ ,4M„ +2) using formula [fTTl fT8l 



dc - 2 



E[e s^ + 2c b)]= «<< r -^^ + 2ch) (14) 

J-°°\ / 2na 2 \/l-2a 2 sl 
the MGF of the filtered current, denoted as *P fs (s) here, can be written as 

2ct 2 s 2 S?(» s ) 

= J] £ 7^ (A=2a 2 X,-) (15) 

f=i (i-w* 

where £>,(f s ) is the ;th component of \b(t s )) in ( fT3l . In this work, we take § = 1/2 for polarized 
noise. The BER can be obtained from ( TT3T > with ifTTl 



±1 r 

2n j Jc± s 



BER >th (/ s ) = — / -AW e -% di) (16) 



where y t h is the detection threshold, + and C+ correspond for 7 SS < y t h, while — and C for 
^ ss > yth- Averaging BERs over all bits in the de Bruijn sequence (N ~ 32 in this work), we 
have (f s = t +sT b , s = 0, ■ ■ ■ ,N- 1) 

N-l 

BER= ^BER yth (f s )/^V. (17) 

s=0 

4. OSNR at the receiver 

For an optical system with ASE power being much larger than other noise sources, the OSNR 
with reference bandwidth B r (O.lnm) can be calculated as 

OSNR .lnm= D P ' r . (18) 
^ASE\tS r ) 

In Eq. dT8b . P s is the time-averaged (noise free) signal power, while PASE(B r ) is the noise power 
within B r . To obtain P s and Pase^B,), one needs to notice that the measurement bandwidth B m 
[e.g., the bandwidth of the transfer function of an optical spectrum analyzer (OS A)] may not be 
the same as B r . Thus the P s in Eq. ( f]~8b becomes the power of the signal filtered by B m , while 
PASE{B r ) becomes the ASE filtered by B,„ and wighted by a factor B r /B m fl9l . 

In a linear optical system, the ASE noise along the fiber can be treated as AWGN (additive 
white Gaussian noise). Thus, for the system of Fig. [4] its OSNR can be simply calculated as 

OSNR Lfi .inm = D P ' ; x ?2 (p AS E(B m ) = [GN in +N (K+ l)]B m ) , (19) 

"ASE{B m ) B, V / 

where PASE{B m ) is the ASE power within B„, and A^o is given by Eq. (1281 . In Eq.([T9l, the filter 
(B m ) effect on the ASE has been neglected. 



CO 



10" 



10" 



10" 



10 



io- 



" - 




CMM 

RK4fl* - ^ - 








- 














































13 15 

OSNR @ O.lnm(dB) 



17 



Fig. 1. BER versus received OSNR for a 20Gb/s 20-span RZ-DPSK system with <J>ai = 
0.2k. Solid: obtained using CMM of Ref. |T). Dashed (dotted): improved CW approach 
of Ref. [5| with CW power being peak power (average power), respectively. Dash-dotted: 
RK4IP approach. All curves, except the dash-dotted, are obtained from Fig. 8 of Ref. J5J. 



In a nonlinear optical system, the ASE noise "amplified by" PG cannot be treated as a white 
noise. Similar to the noise-noise beating given in Eq. ([T3l l. the measured ASE power will only 
relate with the self -beating terms of the noise, which yields 



OSNR 



Ps 



B„ 



NL ,0. \nm 



(20) 



x [P ASE (B m ) = Tr(Oj n PGO m )a 2 = Tr{PGO m O T m )a 2 

Here O m is the low-pass transfer function of the bandpass filter (bandwidth B m ). In Eq. d20l ), 
a 2 and PG are given by Eq. ( T39b and Eq. (O, respectively. 

In the case of traditional OSA-based out-of-band OSNR monitoring, the ASE power can be 
interpolated using 



P ASE (B mi ±AX) = 



Tr 



PG(0 m {-A\)0 T m {-AX) + 0, n {+Al)0l{+M) 



(21) 



where <9„,(±AA) is the filter function centered at ±AA. When AA = 0, Eq. (l2"TT i returns to the 
ASE power in Eq. (|20]), where P A se(B,„) = P ASE (B m ,0). 

5. Applications to DPSK systems 

To show that the new approach to get the noise propagator is numerically applicable, we will 
compare our RK4IP results with the CMM results given by Ref. [5| and with the experimental 
data given by Ref. [7 1. Both consider systems with R/, =20Gb/s, using RZ-50% DPSK modula- 
tion. In the receiver, the optical filter is Gaussian type, while the electric filter is the fifth-order 
Bessel type. In the following calculations, we set 7b = NT), by changing /i in Eq. ( fTTT i. 

5.1. Comparison with CMM results 

We consider the 20-span DPSK system discussed in Ref. @, where BERs using the CMM 
and the (improved) CW approaches were plotted against the received OSNR in the Fig. 8 of 
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Fig. 2. BER versus received OSNR for the multi-span RZ-50% DPSK systems with <S>n = 
0.9. According to Fig. 7 (b) in Ref. (7), where the experimental curves for 5-, 10-, and 
25-span systems were almost the same, here we replot them using a thick solid curve. 



Ref. 0. In fact this system is basically the same as the one shown in Fig.|4]in the Appendix, 
provided that one removes the pre- and postcompensating fibers and their amplifiers in the Fig. 
2 of Ref. and removes the first EDFA and Nj„ in our Fig. [4] Thus the first term of PG [in Eq. 
d2J] needs to be ignored. Like 01201, where OSNR was calculated in the absence of PG, we 
obtain OSNR from Eq. (fT9l l with Nj„ = and (K + 1 ) being replaced by K. According to Ref. 
0, we change OSNR by changing the n sp in Eq. (f28t . As plotted in Fig.|U each span contains 
a transmission fiber followed by a dispersion-compensating fiber (DCF). The transmission fiber 
is / =100 km long with its CD parameter D tx = 8 ps/nm/km. Each span is fully compensated. 
The nonlinear phase accumulated in the fiber, generally defined as = J yP(z)dz with P(z) 
being the time averaged signal power, is 0.27T. The bandwidth of the optical (electrical) filter in 
the receiver is B„ = 1.8-Rj, (B e = 0.65/?/,), respectively. 

To let our results be reproducible, we provide, as detailed as possible, other related param- 
eters below. The DCF in each span is 8 km long with Docf = —100 ps/nm/km. Transmission 
fiber and DCF are assumed to have same fiber loss (a =0.2 dB/km) and same nonlinear coeffi- 
cient ( 7 =2.0 /W/km). The EDFA in each span is used to compensate the total loss in the fiber 
of L = (100 + 8) km. Therefore the signal power at the input of each span, P,„, keeps constant. 
Ignoring the nonlinear phase contribution of DCF Oil , we set P,„=0.7307 mW, obtained from 
4>iV = KYP in {\ - e- a ')/a = 0.2k with K = 20 and / = 100 km. The OSNR is obtained using 
Eq. d} with B m /B r = \35. 

As shown in Fig. Q] the curve using the proposed RK4IP approach (dash-dotted) agrees 
very well with the CMM curve (solid) given by Ref. 0. In Fig. [T] the curves using improved 
CW approach with CW power being transmitted peak power (dashed) and average power 
(dotted) are plotted for comparison. 




5.2. Comparison with experimental data 

The optical system discussed in Ref. can be modelled by Fig. [4] except that each EDFA 
should be replaced by an EDFA followed by an optical filter (Gaussian) with bandwidth of 
5 nm. Also, the input noise Ni n needs to be filtered by an optical filter (9,„ (Gaussian, 5 nm). As 

a result, in Eq. 0), the noise propagator p n (jyk+ \ )L,kLj (k= 0, • • • ,K— 1) should be replaced 

with Oi k p„((k+ \)L,kL\,P n (K,K) = 1 with P„ (K,K) = 0, k , and Gaf n P„{K ,0)P^ {K ,0) with 

Gaf n P n (K, 0)O in Ol n pT(K,Q). Since we only consider the BER vs RX-OSNR curves for the 
multi-span systems plotted in Fig. 7 (b) of Ref. 0, each fiber (L km long) in Fig. [4] contains a 
SMF followed by a DCF In our calculation, all the related fiber parameters are same as those 
given in Table 1 of Ref. 0. In the receiver, the bandwidth of the optical filter is 1 .87/?/,, while 
the bandwidth of the electrical filter is 0.75/?/,. 

We first consider the back-to-back case. Similar to Ref. [7|, we modify the 20Gb/s RZ- 
50% signal at the transmitter by comparing its calculated spectrum l22l and its measured 
spectrum 0]. Due to that the input noise Nj n is filtered by 0,„, the OSNR is calculated using 
Eq. (f20]l with B m /B, = 0.95, yielding the back-to-back RK4IP curve shown in Fig. [2j 

For the 5-, 10-, 25-span systems, their accumulated nonlinear phases are calculated according 
to Eq. (48) of Ref. (0). Because of the spectral modification of the input signal, the optical 
power at the input of each span P( n = Psmf is smaller than Eb/Tb, where Eb is the energy per bit 
before the spectral modification. For example, to get nonlinear phase <t>A? = 0.9 for the 25-span 
system, the fiber input power P;„ = Psmf should be 1.516 mW, which means Eb/Tb =0.127 
mW or GEb/Ti, = 2.316mW (G = 18.197). Different from the DPSK receiver shown in FigH 
where the delay is 7J, = l/Rb =50 ps, the delay in the receiver of Ref. was T fc '=(24.84 



GHz) '=40.26 ps. Thus, the DPSK phase factors given in Eq. ( [Tol l should be modified as 

u Vl — 2 ' m ' m ~ 2 ' ' _ 2 ' 

with N' = N(T b /Tl), T' h = T b +AT b . In Eq. <|22), the phase shift discussed after Eq. ® is 
introduced as 

A G m = 00 + few = ^0 + arctan(VOWfl) . (23) 

To get all the RK4IP curves plotted in Fig. [3] we set 0o = — arctan(vlO) = —1.26, whereas the 
ASE power is calculated using Eq. (|2"TT i with AA = 2B,„. Here, 0o is a calibration constant that 
basically shifts the RK4IP curves in the OSNR direction, while arctan(\/<9SAfR) determines 
the slope of the RK4IP curves. To show this, we plot in Fig. [3] the RK4IP results for the 25- 
span system with Agm = — 1 .26 + arctan( \/OSNR) and Agm = 0. Also, we consider the RK4IP 
curves using Eq. d20l ) to calculate ASE power. Our results for the 5-span, 10-span, and 25-span 
systems confirm that there is almost no difference between the the curve using Eq. (f2TT > with 
00 = —1.26 and the curve using Eq. d20l with 0o = — 100. How to evaluate 0o for a given 
OSNR monitoring technique is expected to be studied elsewhere. 

6. Summary 

For an OFC system with low BER, e.g., < 10~ 9 , efficient way to evaluate the BER is to calculate 
the MGF of the filtered photoelectric current using KLSE. Within linear perturbation, the noise 
impact on MGF (and BER) can be obtained by calculating the noise propagator. 

Different from the CMM approaches in Refs. ifTl fTOlfTTI . a new approach is proposed here 
to get the noise propagator from the accurate LNE using RK4IP method. This approach does 
not need the jitter separation required by CMM. To numerically verify this new approach, we 
consider a 20-span RZ-DPSK system discussed in Ref. |5 |. The BERs obtained using this new 
RK4IP agree well with those using CMM in Ref. 0. 

As discussed after Eq. ([8]), conventional noise propagator obtained by factoring the PG has 
an uncertainty caused by a phase shift in noise. By comparing our results with the experimental 
data of Ref. J7), we find that this deterministic shift is related with the Gordon-Mollenauer 
phase noise variance and it will affect the BER significantly. Further study on such phase shift 
is expected to provide a new approach to reduce the impact of the noise-signal interaction on 
OFC systems. 

Appendix: Accurate LNE in the EDFA-based systems 

The optical field u(z,t) in a fiber satisfies 

Tz =J —M + — W- J ^ U -2 U > (24) 

where a is the fiber loss and Pa,® — d 2 fi /da 2 relates to the CD parameter D(X) (ps/nm/km) 
with p ma ^¥g&l (c=3xl0 8 m/s). The slope parameter p a(om = (^-) 2 [2XD(X) + X 2 D' (X)} 
[D'(X) = dD j^ (ps/nm 2 -km)] can be neglected if bit rate Ri, satisfies Ry > IPmm/Pmwal El- 
Introducing the transformation u{z,t) = v(z,t)e~ az l 2 , Eq. d24l) can be reduced as 

Tz= J — ^ + —1^-^ ™ v - (25) 

In a A'-span system amplified by (K+ 1) EDFAs (cf. Fig.lUi, Eq. $25[ can be modified as 
dv .j5 mm d 2 v Pmmio <? 3 v . _„ 2 

T z =J ^Tdfi + — dfi- je ^Hv-Mz, f) , (26) 




Fig. 4. Low-pass equivalent optical model. The receiver consists of optical and electrical 
filters and DPSK balance detection. Discussions here (except those detailed in Sec.|5j are 
based on the assumptions that iVbfc = No and = G [k = 1,2, ...(K+ 1)]. Also, the fiber in 
each span is assumed to have same length (L km) and same loss (a dB/km). 



where w(z,t) is the ASE forcing modeled as the complex AWGN with correlation 

K+l 

r(z^,t,t')=E{w(z > t)w*(z!,t')} = 8(t-t')8(z-z') 1 £Nok8(z-(k-l)L). (27) 

k=l 

In Eq. d27K the fiber length in each span is assumed to be L (km) long, According to Wiener- 
Kinchine theorem 11231 . Nok in Eq. (|27T > is the ASE PSD (in one polarization direction) at the 
output of the kth EDFA. Suppose each EDFA has the same gain G and spontaneous-emission 
parameter n sp , we have ifTTl 



N k=No=n sp (G-l)h(0 [k= l,2,-K,(K+l)J (28) 

Decomposing optical field v(z,t ) in the fiber into noise-free field Vo(z,t) and its perturbation 
8v(z,t) [i.e., v(z,t) = vo(z,f) + 8v(z,t)] and assuming that |vo| >> |5v| (so that the nonlinear 
terms of 8v can be neglected), Eq. d26l ) can be decomposed as HI 

dv .Pa>a>d 2 Vo Pa>wa> d 3 v Q . 2 

^ = J —^F + —^- je 7|vo|vo (29) 
d8^ =J p^3^ + ^^_ i2e -^ Vo|25v _ ; . e -^ 5v *_ >(Zj/) . (30) 

Noise equation Eq. ( f30b differs from common equations in that the real and imaginary parts of 
the complex noise field need to be treated separately (TJ. 

Denoting a\ = a((0/) = / 8ve~i m ''dt and [M v ] lm = V/_„„ [M^] h „ = p: l+m with 

v, = v(coi) = e- az J \v \ 2 e- jm i'dt, pL, = ;u(cO/) = e~ az / vle'^dt, (31) 
in frequency domain, Eq. ( f30b has the form 

^0/ .Pram 2 .Pmcow 3 , w 1# i . / Nr., i * 

-j- = -j—cofai- j—j—cofat- j2y(z)[M v }i m a m - jY(z)[Mp]i m a* m - jW h (32) 



where Wi = W(z, ft)/) is the Fourier component of the forcing term w(z,t) in Eq. (I26l l. As indi- 
cated in Ref. JTJ, since |vo| 2 in d3~TT l is real, V/ = V*/. SoM v in d32l is Hermitian, or, its real part 
My is symmetric, while its imaginary part M' v is anti-symmetric. Also, as [M^,,, = jU*+ m , both 
the real (M*) and imaginary parts (M^) of are symmetric. 
The matrix form of Eq. d32t is 

da - 

— = La + va + pa* - jW (33) 
dz 

L = jL CD , v = -2jr(z)(M^+jMi), ^ = - J y(z)(M^ + jM I ll ) (34) 

Introducing a = (aR,aj) T (for a = aa + jaj) and W = (Wjr, — Wr) t (for — jW =Wi — ]Wr), Eq. 
(l33l is equivalent to 

fl =(L + v + A)5 + W (35) 



L C d 





vaa 


-v 55 \ 


. A= | 








, v 5S 


vaa y 




-MR y 



^ L CD 

with (IcD)y = -[%i» 2 + %^co 3 ]5 ;> = 2yM£, v ss = -2yM«, HR = 7^, and & = 
— yM^. According to the discussion given below Eq. d32b , V (£l) in Eq. (f36t is antisymmetric 
(symmetric), respectively. Calculation of the Kerr term (v + /i) according to Eq. (l36l l has the 
computational complexity much less than 0(N^), where Nw is the number of Fourier com- 
ponents used for signal representation. In fact, the computational cost of this way is basically 
determined by the FFTs in Eq. d3~Tl i. which has the computational complexity of 0(NwlogNw)- 
In frequency domain, ASE correlation relation (|27| | has its matrix form ( W(cOt) — > W/y/Af ) 

£{W,(z)w;(z')} = S.. : S Lr £ ^r8 z>kL (1 = 1,- ,4M„ + 2), (37) 

where Eq. (1281 has been used. In Eq. (l37l i. 7b = 1/A/ and M n are given by Eq. (fTTl i. Eq. (l37l i 
means that Eq. (l35T l can be equivalently replaced by 

^- = {L + N)a; (N = V + fi) (38) 
dz 

with boundary condition l24l 

£{W(z/)W*(z/)} = ^7=(7 2 7 (39) 

with 7 being a (4M„ + 2) x (4M„ + 2) unit matrix and a 2 being the variance of the real or 
imaginary part of input ASE. 
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